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Abstract 

A coupled channel model is used to study the nature of the scalar mesons 
produced in the decay <f> — > ^tttt. The KK molecular picture of /o(980) is 
found to be in a good agreement with the recent experimental data from 
SND and CMD-2. The structure of the light scalar mesons is elucidated by 
investigating the S-matrix poles and the qq spectral density. 



1 Introduction 

The production of scalar mesons in radiative decays is a valuable source of infor- 
mation on hadron spectroscopy. It was argued jl],j2| that the branching ratio for 
— > 7/0 (980) can be used to make a unique choice among different models of 
/o(980): a conventional quark-antiquark state, an exotic qqqq state ||, and a KK 
molecule ||. Such a possibility to resolve the long debated problem of the / (980) 
structure using a single partial width looks very attractive. However, it was argued 
in m that the dependence of the theoretical predictions for the — > 7/0 (980) decay 
on the /o(980) structure is partly due to differences in modeling. Recent measure- 
ments of the radiative decays — > 77r°7r° by SND || and — > 77r + 7r~~, 77r°7r° by 
CMD-2 ]7[p in Novosibirsk have made it possible to confront alternative models of 
the light scalar mesons with experimental data. 

The goal of this paper is to reanalyze the calculation of the decay width for — > 
7/o(980) in coupled channel models where the /o(980) state arises as a dynamical 
state (a KK molecule which may also have a substantial admixture of a quark- 
antiquark component). Since the meson is nearly a pure ss state, the decay 
— > 77T7T is an OZI-rule violating process which is expected to proceed via a two- 
step mechanism with intermediate KK states. Therefore this decay is well suited 
for probing the KK content of the scalar mesons. The KK molecular state was 



originally proposed in the potential quark model P,|9T,|10||. A dynamical state close 
to A"A^ threshold is also introduced in the coupled channel models [fLl]-|14| and in 
the meson exchange interaction models [13, H3|. A state strongly coupled to the 



ss and A"A" channels near the KK threshold was as well found in the unitarized 
quark model [I7[,[I5]. The coupled-channel model derived from the lowest order 
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chiral Lagrangian |TJ5|,|2t| produces a scalar state dominated by the KK channel. 
General discussions of the nature of the /o(980) can be found in |10|, |13|, |17l, pTh|23 
and references therein. 

Our approach is based on a coupled channel model (CCM) for the tttt and KK 
systems which is similar to the one studied in [13]. The calculation of the decay 



— > 7/0 (980) for point-like particles is summarized in Sec. [2| The details of the 
coupled channel model are given in Sec. |3], and the model parameters are determined 
from a fit to the tth scattering data. The reaction — > 77T7T in a CCM framework is 
studied in Sec. £|. The analytic structure of the Tin and KK scattering amplitudes 
is investigated in Sec. ||. The mixing between the two-meson and qq channels is 
discussed in Sec. ||. In Sec. [7| the physical properties of the scalar mesons in the 
model proposed are discussed and compared to other approaches in the literature. 
The details of the formalism are collected in the Appendices. 



2 The Decay <\> 7/0 

For the benefit of the reader, we begin with a brief summary of the formulas de- 
scribing radiative transitions between vector and scalar states. The amplitude of 
the radiative decay into the scalar meson fo has the following structure which is 
imposed by gauge invariance: 

M(0- 7 /o) = 4e;(p v q,-g u ,(p-q))H(p 2 ,(p-q) 2 ) (1) 

where (e^, p) and (e 7 , q) are the polarizations and four-momenta of the and 7, 
correspondingly, and H(p 2 , (p — q) 2 ) is the scalar invariant amplitude which depends 
on the invariant masses of the initial and final mesons. The polarization vectors 
satisfy the constraints • p = and e 7 • q = 0. Using the three-dimensional gauge 
e 7 = (0, Ery) in the center-of-mass system (CMS) one gets the amplitude 

M(0^ T /o) = {e^) m<l> uH{p 2 ,{p-q) 2 ) (2) 

where is the mass, lj is the photon energy in the CMS, and and e 7 are the 
and 7 three-dimensional polarization vectors in the CMS, respectively. 

For point-like particles, the amplitude M(0 — > 7/0) is given by the diagrams in 
Fig.|TJ (see H and references therein). While the diagram Fig.|I](a) corresponding to 
loop radiation is overall logarithmically divergent, its contribution to the p^q^ term 
in Eq.(|l|) is finite. Gauge invariance enforces the appearance of the seagull diagram 
Fig.|I](b) which contributes only to the g Vi Xp-q) term in Eq.(|l]). Since the sum of the 
loop-radiation and seagull terms is gauge invariant, one obtains the total amplitude 
M(0 — > , -fS) by calculating only the p v q^ term of the loop radiation. The result for 
the scalar invariant amplitude is [^] 

H{p\{p-qf) = e -f^I(a,b) (3) 

a-'4 • (4) 
m z K m z K 
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Figure 1: The minimal gauge-invariant set of diagrams for the decay M(<p — > 7/0) 
in the case of point-like particles: (a) the loop radiation; (b) the contact term. 



where g$ and Q^kk are the <pK + K and f K + K coupling constants and the 
function I (a, b) is defined in [Appendix A| . The radiative decay width is 

The generalization of Eqs. (P],B10) to the case — > 'jtttt where the tttt system has 



total angular momentum J = and isospin I = is straightforward [p5[ : 

M (0^ 7 7T7r) = e^ip^-g^ip-q^H^ip^ip-q) 2 ) . (7) 

Here the scalar invariant amplitude H^-nip 2 , (p — q) 2 ) is given by (compare with 
Eq.©) 



Hnn{p , (P ~ q) ) = 9 2 2 6) t K+K -- nw (8) 

and tK+K--mr is the J = part of the T-matrix for the K + K~ — > tttt scattering. 
The tttt invariant mass distribution has the form 



dT dT ag$ 



«UW - 2 lM^ = 18#^"° " W"- «IX*-„IX, (9) 

where t° K ^_^ = v6 -,r+7r- is the isoscalar — > 7T7T amplitude and k nn is 

the relative momentum of the pions in the final state: 



MlJA-ml , M^ = (p- q y . (10) 

Equation (||D leads to Eq.(||) in the Breit-Wigner (BW) approximation for the 
K + K~ — > 7T7T scattering amplitude 

f _ 9foKKgfoTrn /-, n 

- M^-(M /0 -zr /o /2)2 1 j 

3 



under the assumption that the integral over the tttt mass spectrum is saturated by 
the narrow resonance with mass M/ and width T f (see |Appendix A| ) . 



According to Eq.(|8]), the hadronic part of the amplitude, which contains the 
information about the scalar mesons, is factored out in the form of the T-matrix for 
the K + K~ — > 7T7T scattering, with both the physical region (M wn > 2m K ) and the 
unphysical region (M n7T < 2m^) being relevant to the — > ^tttt decay. It is known 



from the studies of scalar mesons, see [13, E2, E3L G6f and references therein, that 



the analytical structure of the scalar-isoscalar amplitudes near the KK threshold 
is far from being a trivial BW resonance. Therefore a coupled channel model of the 
K + K~ — > 7C7T scattering is required to describe the decay <p — > ^tttt beyond the BW 
approximation. 

3 The tttt - KK Coupled Channel Model 

To describe the interaction in the tttt — KK system with total angular momentum 
J = ++ and isospin I G = + we exploit a coupled channel model similar to that 



of [fL3]| . The two scattering channels, 1 and 2, correspond to the nn and KK systems 
and channel 3 contains a single qq bound state. The T-matrix, as a function of the 
invariant mass squared s, is defined by the Lippmann-Schwinger equation 

T(s) = V + VG°(s)T (12) 

where G°(s) is the free Green function. 

The interaction potentials are taken in separable form: 

/t;u(s)|1)(1| «i2(s)|l)<2| gi 3 \l)(qq\\ 
V = uia(s)|2)(l| v 22 (s)\2)(2\ g 23 \2)(qq\ (13) 
V <7i3|?g)<l| 923\qq)(2\ / 

where the form factors in channel 1 and 2 depend on the corresponding relative 
three-momentum k: 

(k\l) = aW^^Af (M) 
(*|2> = ^ k ) = k^ ^ 



In [0, the potentials Vij(s) were assumed to be energy independent, and the chiral 
symmetry constraints on the scattering amplitude were strictly imposed only in 
the 7T7T channel by adjusting the strength of vn so that the Adler zero is at the 
correct position. In the present case, it is essential to ensure a correct behaviour of 
the KK — tttt scattering amplitude not only in the physical scattering region but 
also down to the tttt threshold. We found it easier to impose the chiral symmetry 
constraints by using energy dependent potentials^. The energy dependence is taken 

lr rhe unitarization of the lowest order in chiral perturbation theory goes along a similar way, 
see [[L9| and references therein. 
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in the form 



vu(s) = b-g n s (16) 
Vi 2 {s) = -gi 2 s (17) 
V22(s) = -g 2 2S . (18) 



With our choice of interaction ( |13| - |T8| ) the analytical solution for the T-matrix can 
be easily obtained. Further details of the model are given in [Appendix B| . As in |13| 



we assume that the diagonal interaction in the KK channel produces a weakly 
bound state in the absence of coupling to the other channels, thus simulating a 
"molecular" origin of the /o(980) resonance. The state \qq) in channel 3 has a bare 
mass M r > 2m k- 

(a) (b) 




Figure 2: The scattering phase 5® and the inelasticity parameter 7]® for the S'-wave 
tttt scattering vs. the tttt invariant mass. The experimental points are from f27j-p9 



The curves show the fits in our model: fit 1 (the dashed lines) involves only the tttt 
scattering data, fit 4 (the solid lines) includes in addition the data for <p — > , jtt 7t 



The model parameters have been determined from the fit of the tttt scattering 
phase Sq and the inelasticity parameter 77$ in the mass range M n7T < 1.5 GeV, see 
Table [l] and Figs.[|(a,b). Since the fit was found to be only weakly sensitive to the 
form-factor parameter A2, its value was fixed, and only the coupling constants gij, 
the bare mass of qq, and Ai were treated as free parameters. The parameter b in the 
diagonal tttt potential Eq. (|T6|) was used for fine tuning of the tttt scattering length 
which was fixed at a[j = 0.22m" 1 . The best fit of only the tttt scattering data (fit 
1) does not lead automatically to a very good description of the tttt invariant mass 
distribution in the decay — > ^tttt. However, the shape of the — > ^tttt improves if 
the dT/dMnn data are added to the fit as shown in Fig.|| (the details are discussed 
in Sec.|j). As a result, our model provides a good description of the whole data set. 



5 



Table 1: The model parameters obtained from the data fits. The fit 1 is based on 
using only the tttt scattering data (the phase 8® and the inelasticity t]q). The fits 2-4 
include in addition the mass distribution dT/dM n7T in the decay — > 77r°7r°: data 
from SND § in fit 2, from CMD-2 @ in fit 3, and the combined data @,§ in fit 4. 



The value of A2 was fixed in all fits. 



fit 


b 


0n 


022 


012 


Ai 


A 2 


013 


023 


M r 




GeV- 1 


GeV- 3 


GeV" 3 


GeV- 3 


GeV 


GeV 


GeV 1 / 2 


GeV 1 / 2 


GeV 


(1) 


5.251 


3.564 


5.326 


3.741 


0.477 


0.7 


2.670 


0.464 


1.136 


(2) 


6.075 


2.953 


3.945 


-0.121 


0.499 


0.7 


2.708 


0.716 


1.094 


(3) 


5.188 


2.944 


3.683 


-0.270 


0.538 


0.7 


2.576 


0.690 


1.111 


(4) 


5.299 


2.954 


3.725 


-0.341 


0.529 


0.7 


2.588 


0.702 


1.105 



The energy dependence of the T-matrix element tK+K^mr calculated in our cou- 
pled channel model is shown in Fig.|j. In addition to a narrow peak due to the /o(980) 
state, this matrix element has a significant contribution from the lower mass region 
corresponding to the a meson. 

4 The Decay <\> — > 77777 in the Coupled Channel 
Model 

The formulas given in Sec. |2| are valid for point-like particles. When a form factor 
in the KK — tttc vertex is included, the — ► 77r7r amplitude contains an additional 
term arising from the minimal substitution —>■ (k^ — eA^) in the momentum 
dependence of the form factors (see || for details). With an appropriate choice of 
the form factors, the sum of three diagrams shown in Fig. |5| becomes explicitly finite. 
In order to study the influence of the form factors on the results for the <fi — > ^nn 
we use a nonrelativistic approximation for the K + and K~ , which is justified by the 
fact that the most interesting region corresponding to the /o(980) resonance is very 
close to the KK threshold. The electric dipole matrix element is factorized into two 
parts describing the K + K~ loop radiation with gauge invariant complement and the 
final state rescattering K + K~ — > titt, correspondingly. The total scalar invariant 
amplitude has the form 

H$(p 2 ,q 2 ) = Tr^-JxiM^tK+K^M^) (19) 
2ir z m K 

which is similar to the relativistic point-like case defined by Eq.(H) where the function 
I(a, b) is replaced by J\(M n7T ). For the definition of ^(M^) and further technical 
details we refer to |Appendix Q . The parameter A refers to the form-factor de- 
pendence on the relative K + K~ momentum in the K + K~ — > tttt vertex given by 
Eq.fll^). Figure ^| shows the dependence of the electric dipole matrix element on the 
7T7T invariant mass for different values of the form-factor parameter A in comparison 
with the relativistic point-like case. 
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Figure 3: The tut invariant mass distribution in the decay <fi — > 77r°7r°. The curves 
are the results from our model: the dashed line corresponds to fit 1, the solid line 
to fit 4. The experimental points are from M and M. 



With our choice of the form factor (^) there is no substantial suppression of the 
nonrelativistic result^ in comparison to the point-like case for the relevant range of 
A = 0.6 — 0.8 GeV (which corresponds to the data fit in our coupled channel model). 
For larger values of A, the full relativistic treatment is needed as the real part of 
the nonrelativistic result becomes sensitive to the short distance contribution (see 
the discussion in J5|). The dependence of the imaginary part of M(<ft — > 77T7t) on 
the form factor in the mass region close to the KK threshold is rather weak. As a 
result, we can neglect the form-factor dependence in the loop calculations and just 
use the relativistic point-like result given by Eq.(||). It is interesting to note that 
the imaginary part can be obtained from the divergent triangle diagram Fig.||(a) 
alone by using the Siegert theorem which ensures a correct behaviour of the electric 
dipole matrix element. 



5 The Poles of the S- Matrix 

The poles of the S"-matrix in the complex s-plane corresponding to the fits in SectfJ 
are shown in Table [| and Fig.[7]. There are five poles related to the resonances in 
our model. The pole Ma located on the sheet II (Im k\ < 0, lm k2 > 0) is very 
close to the KK threshold. The pair of poles, Mb on sheet II and Mb on the sheet 

2 The stronger dependence of the total amplitude on the form factor found in J^] results from 
the use of the dipole form factor which falls off faster than the monopole form factor used in our 
case. 
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Figure 4: The T-matrix for the KK — nn scattering ( J = I = 0) vs. the invariant 
mass of the tttt system M nn . The model parameters correspond to the fit 4 in Table [TJ. 




Figure 5: The minimal set of diagrams giving a gauge-invariant amplitude M(0 — > 
77T7r) in the coupled channel model: (a) the loop radiation; (b) the contact (seagull) 
term, (c) the form factor term. In the nonrelativistic approximation, the particles 
marked with x are on mass shell. 



Ill (1m k\ < 0, Im k 2 < 0), corresponds to a broad structure associated with the a 
meson. The other pair of poles, Mq on sheet II and Me on sheet III, corresponds to a 
broad resonance above the KK threshold, which can be associated with the /o(1370) 
state (for a more realistic description of the S-matrix above 1.3 GeV additional 
poles are needed which are not included in the present model). The position of 
the resonance poles depends on the coupling constants, and this distinguishes them 
from the fixed poles originating from the singularities of the form factors (|15|). The 
latter are located at k\ = HX\ and k 2 = HX2, their distance to the physical region 
being determined by the range of the interaction. In our model these fixed poles 
approximate the potential singularities which correspond to the left hand cut in a 
more general case. 

The origin and the nature of the resonance poles found in our model can be 
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Figure 6: (a) The function I\ = (a — b)I(a,b), a = m 2 ,/m 2 Kl b = Ml n /m 2 K 
vs. the invariant mass of the tttt system M n7T . (b) Its nonrelativistic analog 
J x (M nn ) ■ {m\-Ml^)jm\ for A = 0.5 GeV (solid lines), A = 0.6 GeV (dashed 
lines), and A = 0.8 GeV (dotted lines). 

Table 2: The resonance poles of the ^-matrix in the complex mass plane M = \fs 



(GeV). The definition of the fits is given in Table |1]. 



Pole 


Sheet 


fit 1 


fit 2 


fit 3 


fit 4 


M A 


II 


0.991 


- i0.024 


0.974 


- i0.020 


0.978 


-iO.015 


0.975 


-i0.017 


M B 


II 


0.455 


-i0.237 


0.456 


- i0.232 


0.469 


- i0.234 


0.465 


- i0.236 


M c 


II 


1.214 


- i0.262 


1.401 


- i0.266 


1.424 


-iO.261 


1.417 


- i0.263 


M D 


III 


0.789 


- i0.487 


0.677 


- i0.234 


0.670 


-zO.266 


0.666 


-i0.227 


M E 


III 


1.403 


- i0.320 


1.384 


-i0.261 


1.408 


-i0.250 


1.400 


- i0.249 



elucidated by studying how these poles move in the complex s-plane when the model 
parameters are varied between the physical case determined by the fit and the limit 
of vanishing couplings in the irn channel and between the tttt and the other channels 
(KK and qq): 

vn(s) -»■ x- Vll {s) , 0<a;<l (20) 
v 12 (s) -> x-v 12 (s) (21) 
9i3 - xV 2 - gvi (22) 
923 - x^ 2 -g 13 (23) 

The diagonal interaction in the KK channel with the physical strength of the cou- 
pling g 2 2 produces a bound state close to the KK threshold with mass m K ^ = 
0.97 GeV. 

Our coupled channel model has only one pole Ma = 0.975 — iO. 017 GeV near the 
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-0.04 L 



Figure 7: (a) The trajectories of the S'-matrix poles (fit 4) in the complex mass 
plane (GeV) on the sheets II and III for the KK — tttt and qq — tttt couplings 
increasing from x = to the physical values x — 1. (b) The magnified region of 
the / (890) resonance, (c) The a meson trajectory in the complex plane of the tttt 
relative momentum k\ (GeV) with (solid line) and without (dashed line) coupling 
to the qq state. 
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Figure 8: The poles of the 5"- matrix corresponding to the /o(980) in different models 

[T^,[3|-|37|]: (a) sheet II pole (Im k x > 0, Im k 2 < 0) and 



(b) sheet III pole (Im k\ < 0, Im k 2 < 0). 



KK threshold, which is sufficient for a good description of the tctt scattering data. 
This pole is directly related to a molecular KK state in the absence of coupling 
to the 7T7T channel. The number of the S"-matrix poles near the KK threshold has 
been discussed in the literature (see fl3|, |22|, |2(| and references therein) for a long 
time. While the relation of the /o(980) meson to at least one S- matrix pole close 
to KK threshold is well established, the exact location and even the number of the 
relevant S-matrix poles is model dependent as demonstrated in Fig.|8|. The models 
based on dynamical input (coupled channel, potential, unitarized chiral perturbation 
theory) [|1^,[12|,|1^,|T^,[T^,|T9|,F3^] usually produce only one pole which can be traced 
to a weakly bound KK state in an appropriate limit of the channels coupling. One 
exception to this observation is a coupled channel model |T!| where the second pole 
near the KK threshold results from the interplay of a nearby qq pole with dynamical 
singularities. We were not able to find a good fit with the same feature in our model; 
the main reason appears to be due to the using of the energy dependent potentials 
in Eq. (|l3|) contrary to the case in [|13|]. However, it is not excluded that two-pole 
solutions can be found with some other parametrization of interactions in CCM. It 
remains to be investigated whether two-pole solutions are compatible with the data 
on the — > / y7T7r decay 

l6| routinely produce the second pole 



The i^-matrix parameterizations p^ , P6| , pl| - 
on sheet III. There is, however, a much larger spread in this pole location than on 
sheet II. As shown above, the reaction <p — > 77T7t allows one to probe the scattering 
KK — > 7T7T both above and below the KK threshold, and the ttti mass distribution 
in <fi — > 77T7T is very sensitive to the 5-matrix poles related to the /o(980). Therefore 
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more detailed experimental data on — > ^yrrir would be very useful for reducing the 
present uncertainty about the analytical structure of the S- matrix in the /o(980) 
region. 



6 The Mixing between the qq and Mesonic Chan- 
nels 




Figure 9: (a) The spectral density of the qq state and (b) the corresponding Laplace 
transforms. 

The mixing of the quark-antiquark states with the open meson channels can be 
studied in the CCM by using the probability sum rule for a resonance embedded 
into a continuum |38[] which is described in |Appendix B| . The spectral function p(s) 



which determines the probability density for the qq component in the scattering 
states, as defined by Eq.(^), is shown in Fig.^|(a). In the case of weak coupling, 
the probability density would be well localized near the position of the bare qq state 
at M r = 1.1 GeV. For the physical case, we find a broad peak centered above the 
KK threshold in the region corresponding to the /q(1370) resonance (the poles Me 
and Mc)- There is also a sizable contribution to the qq spectral density from the 
low-mass region of the a meson. The p(s) distribution in the /o(980) resonance 
has a characteristic dip-bump structure resulting from an interplay of the resonance 
pole Ma with a nearby zero of the mass operator II(s). The position and the width 
of the p(s) distribution indicates that an essential contribution to the saturation of 
the sum rule (|52|) comes from the pole Mb related to the a meson, while the narrow 
structure associated with the pole Ma alone plays a minor role. The fact that qq 
coupling with the irn channel significantly enhances the spectral density p(s) in the 
region of the a meson is related to the interplay of the 5*-matrix poles demonstrated 
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in Fig.[7|c: the pole Mb corresponding to the a meson is pushed towards the tttt 
threshold by the pole Mc originating from the qq state. 

Using the spectral density p(s) we can calculate the contribution of the qq scalar 
mesons to the QCD sum rule related to the scalar quark condensate |39| . Figure [|b 



shows the Laplace transform of the spectral density which is used in the sum rule 
analysis: 

I n (M 2 ) = / s n e~ s/M p{s)ds (24) 



The advantage of CCM with respect to the previous studies (see and references 
therein) where the contribution of the qq mesons was usually approximated by one 
narrow resonance is a more realistic shape of the qq spectral density in the low- 
mass region which is important in the Laplace sum rules. As a result, the momenta 
Iq(M 2 ) and I\(M 2 ) have quite different slopes in their M 2 -dependence as shown in 
Fig.|(b). 



7 Discussion 

The theoretical calculations of the decays — > ^tttt and — > 7/o(980) and the cor- 
responding experimental data are summarized in Table |3|. The predicted branching 
ratios for the 77r°7r° and 77r + 7r + channels should be compared with the experimental 
data with some caution. While these branching ratios can be easily defined theoret- 
ically (r^_ >77r +^.- /T^ry^o = 2 for the isoscalar tttt states, e.g.), the data analysis 
relies on the modeling of competing reaction mechanisms. In particular, the interpre- 
tation of the experimental results for the 77r + 7r + channel involves the consideration 
of the interference between the e + e~ — > <p — > ^tttt and e + e~ — > 7p — > '-/tt-k mecha- 
nisms. In this respect, we wish to emphasize again the importance of using a realistic 
K + K~ — > tttt amplitude instead of a simple superposition of BW resonances. The 
channel 77r°7r° which is free from the p meson contribution in the tttt channel appears 
to be better suited for the study of the tttt invariant mass distribution, although this 
case has some background due to the mechanism — > 7r°p° — > 7r°7r°7. 

Our results are in good agreement with recent calculations within the chiral 
unitary approach both for the total width and for the tttt mass distribution. This 
is not surprising because in both cases the /o(980) resonance is produced mainly by 
the attractive interaction in the KK channel. Our result for the total width r^^,,- 
is close to the earlier calculations of the two-step mechanism with the intermediate 
K + K~ state Hf24 which used the BW approximation (Jg) with the coupling constant 



g 2 QKR /A-K = 0.6 GeV 2 . In our model we can define an effective coupling constant 
9f KK by approximating the KK — > tttt amplitude (f|S|) by a BW resonance which 
leads to slightly higher value of g 2 oKR /4:7i = 1.1 GeV 2 . 

The effect of the form factor in the KKmi vertex was studied earlier in || 
where a suppression by a factor of about 5 was given as an estimate. This signif- 
icant suppression resulted from the use of a very soft dipole form factor with the 
characteristic parameter p = 0.14 GeV (Eq.(4.14) in [p[) which was suggested by 
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Reaction 


Branching ratio 


Ref. 


Theory 


-> 7 /o(980) -> 7^ 


2.6 ■ 10" 4 


E 


v ^ 7,M y ° u J 


10.1 n~ 4 
i.y • ru 


[ 


m 


-> 7/o (980) (point-like) 

— > 7/0 (980) (with form factor) 


1.4- 10" 4 
0.3 • 10" 4 




§ 


— ► 7(/o ± ® ) — * 7 7r7r {q^Q^ model) 
— * 7(/o ± cr) — > 77Tvr (ffi model) 
— *■ 7(/o ± cr) — > 77T7T (ss model) 


~ ru 

~ 10^ 5 
~ 5 • 10" 5 


[ 


m 

12f 
42j 


- 7(vr°vr°) J=0 

-> 7(7r + 7T-)j =0 


0.8- 10" 4 
1.6 • 10~ 4 


I 


25|j 


-> 7(7T 7r°)j =0 

-> 7(vr 7r )j =0 , (AW > 0.9 GeV) 
-> 7(7r + 7r-), /= o 


1.2 ■ 10- 4 
0.34- 10~ 4 

2.3 • 10" 4 


this paper 
this paper 
this paper 


Experiment 


- 7(vr°vr ) J=0 

^ 7(vr vr ) J=0 , (M^o > 0.9 GeV) 
- 7 /o(980) 


(1.14 ±0.10 ±0.12) • 10~ 4 
(0.5 ±0.06 ±0.06) • 10" 4 
(3.42 ±0.30 ±0.36) • 10" 4 


i 
1 
1 


- 7 /o(980) 

— > 77r°7T 

-> 7vr vr , (AW > 0.9 GeV) 


(1.96 ±0.46 ±0.50) • 10" 4 
(1.08 ±0.17 ±0.09) • 10~ 4 
(0.57 ±0.06 ±0.04) • 10- 4 


E 
1 
i 



Table 3: The calculated branching ratios for — > 7(7T7r)j =0 and — > 7/ (980) in 
comparison with the experimental data. 



the study of the — > 77 decay |41| . However, the 27 decay is related to the short 
distance behavior of the KK wave function, while we are interested in the <fiKK 
vertex at moderate relative momenta (in general, there is no unique relation be- 
tween these two properties, see the discussion in ]4"I|] ). It is difficult to justify such 
a low value of the form-factor parameter in the CCM: a good fit of the scattering 
data requires much harder form factors as discussed in Sec. |3|. A weak effect of the 
form factor found in our case is consistent with the results in |J if the form-factor 
parameter \i > 0.5 GeV is used there. 

In view of this overall agreement between different calculations of the two-step 
mechanisms with the K + K~ intermediate state we find it difficult to agree with the 
statements P, fi2]| that the KK model of /o(980) can be excluded on the ground 
of its alleged conflict with the — > 7/0 (980) data. Since the hadronic part of 
the amplitude — > 77T7T is factored out as the matrix element t K ^_ ww or, in the 
simplified case, as the coupling constant g faKR > the self-consistency of the KK 
model of / (980) can be examined by analyzing the value of g j kk- F° r this purpose 
we use a well known result from nonrelativistic scattering theory: the residue of the 
scattering matrix at a pole corresponding to a bound state is uniquely related to 
the asymptotic normalization constant of the bound state wave function which for 
a weakly bound state depends, in leading order, only on the binding energy. For 
a weakly bound KK state, the relation between the gj KK an d the binding energy 
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Ef, = 2mK — rrif has the form 

9%kk 



An 



32m K K(l + 0(K/X)) (25) 



K = Jm K E b (26) 



where A is the range of the KK interaction. Taking E b = 4 MeV and neglecting 

a 2 - 

small corrections of the order k/X one gets Js ^ L _= 0.7 GeV 2 which is very close to 
the values discussed above. Therefore a large KK component is naturally expected 
for the /o(980) state regardless of further details of particular models. 



8 Conclusion 

The decay <fi — > '-firir has been studied in an exactly solvable coupled channel model 
containing the hit, KK, and qq channels using separable potentials. The /o(980) 
resonance corresponds to one S'-matrix pole close to the KK threshold; this pole 
has a dynamical origin and represents the molecular-like KK state. The molecular 
picture of the / (980) meson is found to be in a fair agreement with the experimental 
data. We confirm the assessment of U that the earlier conclusions about suppression 
of the cf) —>■ 7/0 (980) branching ratio in the molecular KK model were partly related 
to differences in modeling. 

The lightest scalar meson, a, has a dynamical origin resulting from the attrac- 
tive character of the effective irn interaction, with a partial contribution from the 
coupling via the intermediate scalar qq states. The distinction between genuine qq 
states and dynamical resonances, a and /o(980), can be illuminated by considering 
the limit N c — > 00 where the qq states turn into infinitely narrow resonances while 
the dynamical states disappear altogether. 

The structure of the qq state embedded into the mesonic continuum has been 
analyzed using the calculated qq spectral density. The gross structure of the quark- 
antiquark spectral density p(s) is related to the / (1370) resonance. There is also 
a significant contribution to p(s) in the low mass region (a meson) which is related 
to the strong coupling between the 7T7t and qq channels. The same approach can 
also be used for the QCD sum rules related to the gluon condensate by extending 
the coupled channel model to include the mixing with the scalar glueballs. The 
consideration of this topic is beyond the scope of this paper. 
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Appendix A The function /(a, b) and the decay 

widths 



The function I(a, b) is given by (see e.g. P^J and references therein) 

1 2 + 



/(a, b) 



2(a - 6) (a - 6) 2 
-(arcsin^) 2 



1 l + VIEg -jTr 

4 V 1 - Vl - 4x 2 

\/4a; 2 — 1 arcsin - -t= , 
2v^ 



x > 



±VT^4^ fin i^^S - in) 
2 V 1 - Vl - 4x 2 / 



x > 

X < 



(28) 



(29) 



1 - Vl - 4x 2 

The (j)K + K~ coupling constant is related to the decay width by 



K + K~ 



9 \ ( m * _ A TO 2 N3/2 



(30) 



The relation between the coupling constants and decay widths for the scalar meson 
fo has the form 



r(/o - vrvr) 



1 9%, 



2 167rm 2 Q 



(31) 



where the extra factor \ accounts for the identity of the two pions in the final state. 



Appendix B The Coupled Channel Model 



The free Green function G°(s) is a diagonal matrix: 



G°(s) 





V o 





% 





G° 2 (s) 






G° 3 (s) 



where the single-channel Green functions have the form 



G?(s) 



G° 2 (s) 



G° 3 (s) 



\h)(h 



7r Jo s/4 - (m 2 + k 2 ) 

2 roo 



k\dk\ 



7r 7o s/4 - (m K + k\ 
G° 3 (s)\qq)(qq\ 



(32) 



(33) 
(34) 
(35) 
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Here \k\) and I&2) denote the free tttt and KK states with relative momenta k\ and 
k 2 , respectively. The state \qq) in channel 3 has a bare mass M r . With the form 
factors given by Eq.(JfJ) the matrix elements of the Green functions are 

<5M = WGJMW = 2{k J+a n) ! • » = (36) 

where k„(s) is the relative momentum in the channel n: 



h(s) = y/s/4 - ml (37) 



k 2 (s) = ^s/A-m 2 K . (38) 

The im elastic scattering amplitude f nn (s) and the nn—KK amplitude f^-KR{ s ) 
have the form: 

f m {s) = -(k 1 \T(s)\k 1 ) = -ak 1 ) 2 ^^- (39) 

D(s) 

L,-kk(s) = -(h\T(s)\k 1 ) = -ak 1 )ak2)^- (40) 

D(s) 

where 

D(s) = l-u n (s)G 1 {s)-u 22 {s)G° 2 {s) 

+ (un(s)u 22 (s) - uj 2 (s))G%s)G° 2 (s) (41) 

N n (s) = un(s) - (u n (s)u 22 (s) - u 2 12 (s))G° 3 (s) (42) 

N 12 (s) = u 12 (43) 

un(s) = v n (s) + J ^ (44) 

U12W = + ^ (45) 

U22 ( S ) = V22 ( S ) + _^I|_ . (46) 

The connection between the partial wave S-matrix and the scattering amplitude 
/m is given by 

S^s) = V ° (s)e 2iS o^ = 1 + mMs) (47) 

where 5®(s) is the scattering phase and t]q(s) is the inelasticity parameter. The 
T-matrix element is related to the amplitude (fiO|) by 

^KK-wn = 8iry/sf K K_ W7T . (48) 



The spectral density of the qq state has the form [13 



p(s) = -L(G 3 ( a -i e )-G 3 (s + i e )) , (49) 
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where G 3 (s) is the exact Green function in the \qq) subspace: 

CW = (w\G(sm= s _J_ n(s) (50) 
and II(s) is the mass operator of the qq state: 

m % 2 3 G°( g ) + gl^g) + (2v 12 (s)g 13 g 23 - ^(s) - g^nC*))^ (*)gg(*)) 
U 1 - v 11 (s)G° 1 (s) - v 22 (s)G° 2 (s) + (v n (s)v 22 (s) - i;? 2 (s))G?(s)G§(s) 

(51) 

The spectral density p(s) satisfies the normalization condition 

p(s)ds = 1 . (52) 

Ami 



Equations (^9 , 52) represent the completeness relation projected onto the qq channel 



and the normalization (qq\qq) = 1. 

Appendix C The Nonrelativistic Approximation 

The amplitudes corresponding to the diagrams in Fig. |5| have the form: 

M NR = gg» (/, + h + Q = ^ tK : K ;~7 ^ UM^) (53) 



where 



d 3 k (eg ■ k)(g 7 ■ (k + q/2))F(|k + q/2|) 
2£ fc (jfe 2 - + ie)((Jfe - g) 2 - m 2 K + ie) 



(54) 



±™ K (e -k)(e 7 -k)F(|k + q/2|) ^ 



it 



m^M^J (A -k 2 + ze)(A„-(k + q/2) i 

or. ^f d " k F (l k +q/ 2 l) r «n 

2 ^-^y ^ -<+.)= (56) 

^ ^ k _?W_ (57) 



M_ 7 (A^ - k 2 + ie) 
rf 3 k (^-k)(s 7 -kA)^F(|k|) 

/c " 2 J2E k (k*-m 2 K + ie) " (5§) 

= Icl+/c2 (59) 

= _(fl^ /A (60) 

Afy j (A<£ - k + ze) 

3M 7 (A^-^ + ze) 2 1 J 



Here E\. = J k 2 + m 2 K , A^ = (m^ — 2rriK)mK, A nn = (M nn — 2m^)m^:, the form 
factor -F(|k|) = ^2(|k|) according to Eq.fllBD, the subscript A in J x (M n7T ) refers 
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to the form-factor parameter. In deriving the nonrelativistic approximation (55- 



9p we keep only the positive-energy parts of the kaon propagators and substitute 
d 3 k/Ek — > d 3 k/rriK in the final integral. The total amplitude M vanishes at 
uj = as expected for the electric dipole transition. The diagram Fig.|5](c) contains 
two terms according to Eq. ([59]) which can be combined with the loop radiation term 
(|55D and the contact term (|5j]) in a way that the combinations (I a +I C 2) and (Ib+I c i) 
are explicitly finite and both vanish at uj = |q| — > 0. The integrals (|55| - |6l| ) with our 



choice of the form factor can be straightforwardly calculated in analytical form. 
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